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Abstract 

Statistical theory of the dielectric susceptibility of polar liquid crys- 
tals is proposed. The molecules are calamitic or bent-core but the per- 
manent dipole moment is perpendicular to the molecule long axis. The 
ordering of the phase is described by means of the mean-field theory 
based on the Maier-Saupe approach. The theory is used to calculate 
the temperature dependence of the order parameters and the suscep- 
tibilities. The phase diagram with four phases is obtained: isotropic, 
uniaxial nematic, uniaxial ferroelectric, and biaxial ferroelectric. Four 
critical points are predicted. 

PACS numbers: 61.30.Cz, 77.84.Nh 

1 Introduction 

Liquid crystals are always built from anisotropic molecules [T]. For a long 
time calamitic or discotic molecules were considered as appropriate molecules 
for the formation of liquid crystal phases. But it appeared that bent-core 
molecules can lead to many interesting phenomena related to polarity and 
chirality [2j. Phases and phases transitions that can take place in the bent- 
core systems were studied by Lubensky and Radzihowsky p]. They showed 
that to completely characterize phases a third-rank tensor order parameter 
is necessary in addition to the vector and the nematic (second-rank) tensor 
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order parameters. Twelve different liquid phases were identified and many 
symmerty-allowed transitions among them were analysed in detail. The bent- 
core molecules were studied also by means of computer simulations [4j-[7]. 

We would like to study a simple model that can capture main features 
of bent-core molecules. Let us consider a system of N molecules contained 
in a volume V at temperature T. We assume that the potential energy of 
the interactions ^(-Ri, -R2) depends only on the molecule orientations Ri and 
i?2- The orientation of a molecule is described by the three Euler angles 
R = {(j),6,ip) or by the three orthonormal vectors {l,7fi,n). The potential 
energy of interactions has the form 

V{Ri, R2) =Vo + ViPiih ■ h) + V2P2{ni ■ na), (1) 

where Pj are the Legendre polynomials, the vector n determine the long 
molecule axis. If f 1 = then the Maier-Saupe theory is recovered, where the 
nematic phase is present for f 2 < 0. The Vi term breaks Dooh symmetry and it 
partly describes bent-core molecules with C2V symmetry. On the other hand, 
the f 1 term can be connected with the permanent electric dipol p = pi oi 
a molecule because many bent-core molecules have a transverse dipol which 
bisects the bend angle 0. 

Our aim is to create the theory of the susceptibility which is a tensor in 
the case of anisotropic phases 

- (o^ 

where P and E are the polarization and the electric field, respectively. In 
the first approximation, the polarization can be written as 

Pa = ^{L), (3) 

where p{la) is the average value of the dipole component in the a-direction 
in the presence of the electric field. Two different averages are involved here. 
The bar refers to the potential energy of the dipole moment in the electric 
field and the brackets to the nematic potential. The potential energy of the 
dipole p in the electric field E is p ■ E. The linear approximation is applied 
to resolve the influence of the electric field 

(L) = + PplfsE^)) = f3p{lJp)Ep. (4) 
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Thus the susceptibihty has the form 

^P'^ 11 1 \ K\ 
= y {lalp)- (5) 

The components of the susceptibihty in the nematic phase oriented along the 
z axis are 

Np\l-S) 

60X.. - ^y^^^ , (6) 

Np\2 + S) 

"^"^ WkeT ' ^ > 

where S = {P2{nz)) is the order parameter. More advanced calculations were 
given by Maier and Meier [8] who extended the Onsager theory of the sus- 
ceptibility to nematic liquid crystals. If we neglect the induced polarization 
we can write their results as 

.V - ( ^X + S \ Np\l-S) 

= /3x + 3\ iV(2 + 5) 
'"^"^ [2X + 3) QVkBT ' 

where x is the average susceptibility of the phase. 

The organization of this paper is as follows: In Sec. |2] the mean-field 
theory of the phase ordering is provided. In Sec. [3] the formulae for the 
susceptibility are derived in the case of the nematic and the ferroelectric 
phase. Section H] is devoted to some applications of the presented theory. 
Section [3] contains a summary. Appendix A provides the definitions and 
main properties of the basic functions. 



2 Mean-field theory 

The potential energy of molecular interactions ([1]) can be rewritten in the 
form 

R2) = v^ + v,Ei'^iR,'R,) + V2Ej,fiR,'R,), (10) 

where E^J are the basic functions defined in Appendix A [9j. The basic 
functions will facilitate many future calculations. The molecular orientation 
in the phase is described by the distribution function 

JdRf{R) = l. (11) 
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The mean of any function A — A{R) we calculate as 

{A) = jdRf{R)A{R). (12) 

The state of the system is described by a series of order parameters {E^^}) 
but the most important have j — 1 oy j — 2. The internal energy of the 
system is 

N f 

U ^ - J dRidR2f{Ri)f{R2)V{Ri,R2), (13) 
whereas the entropy of the system has the form 

S = -keN J dRf{R) Hf{R)CN]. (14) 

— * 

The energy of permanent dipole moments in the electric field E is 

Ue^N{p-E). (15) 
The total free energy of the system is the sum 

Ftot = U + Ue-TS. (16) 
In the mean-field approximation a potential energy is 

w{R)=j:j:^i^H'^iR)^ (17) 

j nu 

and we find from the Boltzmann distribution that 

f{R) = exp[-l3W{R)]/Z, (18) 
where Z is a normalization constant. The consistency condition 

W{Ri) = ldR2f{R2)V{Ri,R2) -p{Ri)-E (19) 

leads to equations 

= v,{El!^)5,i-pS,,{EJi^-EyS-,^ + EJo^.): (20) 
<^ = V2{EjS)S.o. (21) 
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It is useful to introduce the dimensionless parameters Sj^J = —(3w^Jj} for j > 



in/(i?) = EE4i^4i^(^) 



(22) 



j flU 



Sqq is responsible for the normalization and it depends on other Sjj'J with 
j>0 



c(o) 

^00 



In 



Vi>o 



Uil^{E^SES)-{E^S){ES), 



dsl'J 



Now the equations (l20l)-( l2T]) have the form 



(3p6ui{EJif, - Eyd^if, + EJo^), 
0. 



Sj^J + (3v2{ElS)S.o 
The solution is orient at ionally stable only if the matrix 



kl 



l>0 P 



(23) 

(24) 
(25) 
(26) 

(27) 



(28) 
(29) 



(30) 



is positive definite. The isotropic phase is orient at ionally stable if /3vi > —3 
and /3t>2 > —5. 

3 Dielectric susceptibility 

The dielectric susceptibility tensor is defined by the Eq. ([2]) and we calculate 
the orientational polarization as 



Pa = {la)Np/V. 



(31) 
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Note that the polarization depends on the electric field via the distribution 
function. The components of the susceptibility are 



^oXxx — 



^oXyy 



Np 



11 "^^Mi , 7- ri2 '^'-'mO 



11 r^i£L I 7-7-1; 



a 0(2) 
12 ^"-"mo 



Np 



(1) 



V^T-rll^^V 



no(2) 
12 '^'^^0 



where the derivatives are calculated from the equations 



55, 



(1) 
Ml 



+ /3t^i 



ao{l) 
rll'^'-'pl 



(2) 



^ 4" ''^ dE. 



(1) 



95 



(1) 
Ml 



/3fi 



(1) 
pi 



no(2) 
22 ^^pO 



dE,. 



„ dE. 



dE. 



21 

HP 



ds: 



(1) 



pi 



dE,, 



Qo(2) 
22 '^^pO 

(2) 



dS 



(2) 



V + V r/12 '^'^po 

MP ^iZT" tip 

dS^^ 



dE, 

d^ 
dE, 



MP ^IZT- ^ /-^ / 



MP 



dE, 



0, 

-(3p8-i^, 
0, 

l^pSoH, 



= 0. 



(32) 
(33) 
(34) 



(35) 
(36) 
(37) 
(38) 
(39) 
(40) 



Now we are in the position to discuss the results for possible phases. 



3.1 The isotropic phase 

For the zero field all order parameters are equal to zero. For the nonzero 
field the phase has the symmetry C^ov (the symmetry of the electric field). 



z 



Np^ 



V{3kBT + vi) ■ 



(41) 



6 



For Vi < we get the Curie- Weiss law describing the divergence of x when 
we approach the Curie temperature from above. For i>i > the susceptibihty 
is finite. 



3.2 The uniaxial nematic phase 

For the zero field the phase has the symmetry and the parameter Sqq 
is nonzero for the phase oriented along the z axis. The nonzero order param- 
eters are (Eqq) for j even. For the parallel field the phase has the symmetry 
Cocv whereas for the perpendicular field the phase has the symmetry 
The expressions C/^^ are nonzero. 

eoXx. - v{l + Pv^Ully ^ ^ 

eoX = ' (44) 

Note that Xxx = Xyy 

3.3 The uniaxial ferroelectric phase 

For the zero field and for the parallel field the phase has the symmetry Coov 
The phase is oriented along the x axis. We have nonzero S^i and S'20'' = 
— a/SS'oo''. The nonzero order parameters are: {E^^), (E^) — —\/3{E^), 
(£"02 ), (-E22 )• III ideally ordered phase the values of the order parameters 
are: {e[\^) = 1, {eS!) = 1/4, {Ej^} = {E^^} = -^/4, = 3/4. The 

most important nonzero elements U^J^ are: UH, Uqq, UI2, UIq, Uq'}, 

tt12 7-7-22 tt22 tt22 tj22 tt22 tj22 
<-^-l,-2) <-^22) <-^20' <-^00' <^11' <^-l,-l) '^-2, -2- 

e,Xxx y y ) {l^pv^Ull)A,-Pv^Ul^A2 + f^v,UliA,' ^ 

A^ = {l + l3v2Ull){l + f3v2Ull)-{l3v2Ull){l3v2Ull), (46) 

A2 = (/3t;2f/2'i)(l + /3^2f/o'o')-(/3^2f/o'i')(/3^2f/2'o'), (47) 
A3 = {l3v2Ull){l3v2Ull)-{(5v2Ull){l + (5v2Ull). 



where 
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3.4 The biaxial ferroelectric phase 



For the zero field and for the parallel field the phase has the symmetry C2. 
The phase is oriented along the x axis. We have independent nonzero S[\ 



(£'02 ), {E20 ), {E22 ), and {E_i_i). The susceptibility is described by Eq. 

iSl). 

4 Results 



In this section we carry out calculations for different physical systems of 
polar molecules described by the considered model. The phase diagram of 
the model is shown in Fig. 1. Four phases are present: isotropic, uniaxial 
nematic, uniaxial ferroelectric, and biaxial ferroelectric. Dashed half-lines in 
the picture denote different physical systems with the fixed parameters Vi 
and V2. Four critical points from Ci to C4 are predicted. Ci = (—3.0, —3.3) 
and C2 = (—3.4, —3.8) are tricritical points where transitions changes from 
second-order to first-order. C3 = (—3.0,-4.1) and C4 = (-2.5,-4.54) are 
points where three phases coexist in equilibrium. 

Let us start from the systems known from the previous studies [lOj. For 
the case of f 2 = [half line (a) in Fig. 1] we get the simple system with the 
second-order transition from the isotropic to the uniaxial ferroelectric phase 
at Tc = —Vi/SkB- For the case of f 1 = [half line (d) in Fig. 1] the interac- 
tions are uniaxial (small dipols) and there is the first-order transition from 
the isotropic to the uniaxial nematic phase. The temperature dependence of 
the inversed susceptibility is presented in Fig. 2. In all pictures, T denotes 
the dimensionless temperature. T = 1 corresponds to the transition from the 
isotropic to the nematic or ferroelectric phase. The susceptibilities are ex- 
pressed in Np'^/{Veo\vi\). On decreasing the temperature the susceptibility 
splits into Xzz < Xxx- Xxx runs to the infinity whereas Xzz remains finite. This 
is typical for the uniaxial molecules with the dipol moment perpendicular to 
the symmetry axis. 

In the system described by half line (b) in Fig. 1 there are the second- 
order transitions from the isotropic phase to the uniaxial ferroelectric phase 
(T = 1) and next to the biaxial ferroelectric phase (T = 0.81). The temper- 
ature dependence of the order parameters and the inversed susceptibility is 
presented in Fig. 3 and 4, respectively. We note that the susceptibility re- 
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veals the transition from the uniaxial ferroelectric to the biaxial ferroelectric 
phase. 

A very interesting situation takes place in the system described by half 
line (c) in Fig. 1. The temperature dependence of the order parameters and 
the inversed susceptibility is presented in Fig. 5 and 6, respectively. On de- 
creasing the temperature we meet the first-order transition from the isotropic 
to the uniaxial nematic phase (T = 1) and the second-order transition to the 
biaxial ferroelectric phase (T = 0.96). In the uniaxial nematic phase we have 

Xxx > Xyy ~ Xzz- 

In the case of the system described by half line (e) in Fig. 1 two phases 
are present: isotropic and uniaxial nematic with the first-order transition 
between them. The order parameter {Eoq ) is the most important and the 
susceptibility is finite. 

5 Summary 

In this paper, we presented the statistical theory of the dielectric susceptibil- 
ity of polar liquid crystals. The molecules were calamitic or bent-core with 
the permanent perpendicular dipole moment. We focused on spatially homo- 
geneous (liquid) phases because the understanding and classifying of those 
phases is a prerequisite to study of more ordered phases, which in addition 
break translational symmetry. 

Liquid phases of bent-core molecules from the considered model are listed 
in Table 1. Phase transition types obtained in the present paper can be 
summarized as follows: 

1. I N transition: first-order; 

2. I V transition: second-order or weakly first-order (second-order in 

my, 

3. I ^ V + 2 transition: first-order (not present in [3j); 

4. N ^ V + 2 transition: second-order; 

5. V ^ V + 2 transition: second-order or weakly first-order (second-order 

in m- 
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Every transition has unique features which can be observed from the temper- 
ature dependence of the order parameters or susceptibihties. Generally, the 
classification by Lubensky and Radzihovsky [3] is completed. The results for 
the transitions I ^ V , I ^ V + 2, and V ^ V + 2 was obtained also by Met- 
tout et al |13| in a theory with two vectors representing the six-dimentional 
order parameter associated with the transition from the isotropic liquid to a 
polar nematic phase. 

Polar nematic phases of thermotropic liquid crystals may have interesting 
technological applications, and that is why a number of theoretical studies 
have focused on the possibility of realizing such phases |T3]-[T^. One may 
expect for such phases, a rich variety of unusual textures and defects. The 
existence of biaxial and uniaxial ferroelectric phases requires further experi- 
mental confirmation. 

An extended version of the considered model can be obtained from the 
Straley model of biaxial nematics |12] by the inclusion of our Pi{li ■ h) term 
or the Pi{ni ■ 722) term. Then the biaxial nematic + 2 phase with 
symmetry would be present and new transition types would be possible: 
+ + oriV + 2-^y + 2. 

Appendix A 

Below we list the properties of the functions -E^i)- The functions can be used 
to describe any physical quantity which depends on the three Euler angles. 

1. The definition is 



where R = (0, 6, ip) (the three Euler angles), j is a non- negative integer, 
fi and z/ are integers. Functions Dj^^ are standard rotation matrix 
elements yuj and 



EjiHl^ = {j^y^''^^''^ + .) + (!- .)sign(^)sign(.)] 



+ sign{u){-iyDlil,{R) + sign(/i)(-l)^DL^i,(i?)], 



(49) 




for X > 
for X < 0. 



(50) 



Note that 



sign(— a; 



sign(x) + 26ox- 



(51) 



10 



2. The functions E^^J are real. 

3. The functions satisfy the orthogonality relations 

/ dRE^}{R)EfJ{R) 

4. Let us assume that the three Euler angles R = (0, 6, ip) determine the 
orientation of the three unit orthogonal vectors (/, m, n). The functions 
E^^} can be expressed in terms of the vector coordinates. 



E«(/2)=n„ 


(53) 


E^^{R)^h, 
E[^{R)^n,, 

EiiiR) = Ix, 


(54) 
(55) 
(56) 


E!S{R)-l{-l + ^nl), 


(57) 


Eg\R)^^{-l + nl + 2ll), 


(58) 


Ei^{R) = ^{-l + nl + 2ni), 


(59) 


Eif{R) = ^(-3 + nl + 2ll + 2nl + All). 


(60) 
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Figure 1: Phase diagram of the model considered in the paper, where Vi and 
V2 are the parameters of molecular interactions, = kBT. Four phases 
are present: isotropic, uniaxial nematic, uniaxial ferroelectric, and biaxial 
ferroelectric. Four critical points from Ci to C4 are predicted. Dashed half- 
lines denote different physical systems: (a) V2 = 0, (b) V2 = Vi, (c) V2 = 2vi, 
(d) Vi — 0, and (e) V2 — —2vi. For a given physical system on decreasing the 
temperature we are moving from the center (0,0) to the edge of the figure. 
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3 - 
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I 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 

T 



Figure 2: Temperature dependence of the inversed susceptibility for vi — 
[half line (d) in Fig. 1]. On decreasing the temperature the susceptibihty 
sphts into Xzz < X XX- 



Table 1: Liquid phases of banana-shaped molecules, their symmetries in 
Schoenflies notation, and nonzero variables in the considered model. Four 
phases are present: isotropic (/), uniaxial nematic (N), uniaxial ferroelectric 
(y), and biaxial ferroelectric {V + 2). 



Phase Symmetry Variables 
— 0(3) T 

V CooV Sli, S2Q — —^/SSq^ 
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0.2 0.4 0.6 0.8 1 1.2 

T 

Figure 3: Temperature dependence of the order parameters for vi = V2 [half 
hne (b) in Fig. 1] . There are the second-order transitions from the isotropic 
phase to the uniaxial ferroelectric phase {T — 1) and next to the biaxial 
ferroelectric phase (T = 0.81). 
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0.6 0.7 0.8 0.9 1 1.1 1.2 

T 



Figure 4: Temperature dependence of the inversed susceptibility for vi — V2 
[half line (b) in Fig. 1]. 
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-0.2 I 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 

T 

Figure 5: Temperature dependence of the order parameters for V2 = 2vi 
[half hne (c) in Fig. 1]. There is the first-order transition from the isotropic 
phase to the uniaxial nematic phase (T = 1) and the second-order transition 
from the uniaxial nematic phase to the biaxial ferroelectric phase (T = 0.96). 
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Figure 6: Temperature dependence of the inversed susceptibility for V2 — 
[half hne (c) in Fig. 1]. 
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